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ABSTRACT 

We study the mass-density relationship n oc m x in molecular cloud condensations 
(clumps), considering various cquipartition relations between their gravitational, ki- 
netic, internal and magnetic energies. Clumps are described statistically, with a density 
distribution that reflects a lognormal probability density function (pdf) in turbulent 
cold interstellar medium. The clump mass-density exponent x derived at different 
scales L varies in most of the cases within the range —2.5 < x < —0.2, with a pro- 
nounced scale dependence and in consistency with observations. When derived from 
the global size-mass relationship m oc Z 7 s lob for set of clumps, generated at all scales, 
the clump mass-density exponent has typical values —3.0 < x(7 g iob) ^ —0.3 that de- 
pend on the forms of energy, included in the equipartition relations and on the velocity 
scaling law whereas the description of clump geometry is important when magnetic 
energy is taken into account. 
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1 INTRODUCTION 

Observational evidence testifies that stars are born in 
molecular clouds (MCs) which sizes and mean den- 
sities vary in large ranges: 0.1 — 100 pc and 1 — 
10 5 cm -3 , respectively. T ypical sites of star formation 



are d ense regions in MCs (Enoch et alj|2007l : lAndri" et al.l 
20101), often assoc i ated with young stellar objects (e.g., 
Johnstone fc Baifvl |2JD06|), wherein local gravitational in- 
stability leads to collapse and/or fragmentation and for- 
mation of protostellar cores. According to the modern 
paradigm in the star formation (SF) theory, such condensa- 
tions r esult from shock s , generated by supersonic turbulent 
flows (lElmegreenlll993l: [Klessen, Heitsch fc Mac Low] l2000l: 



Padoan fc Nordlundl |2002| ; 

2009^ 3 



Kritsuk et ail 120071: iDib et all 

Mac Low fc Klessen 2004 for 



20081 : Kritsuk et al 
a review) . Q 

Two main results f rom the theory of in compressible tur- 
bulence, developed by iKolmogorovl (|l94ll ). are the scaling 
law for the velocity dispersion and the kinetic energy spec- 
trum. Prominent characteristics of turbulence in the cold 
interstellar medium (ISM) are compressibility and super- 
sonic flows which lead to shocks and discontinuous density 
fluctuations. On the other hand, the hierarchy of scales and 
the cascade of the kinetic energy in them, predicted in the 
Kolmogorov's theory, are present in the cold ISM as well, 
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although with different scaling relations. O bservational e vi- 
dence are "Larson's first and second laws" |Larsonlll98ll ): 



u oc L 13 
(n) oc L° 



(1) 
(2) 



where L is the spatial scale, u is the rms velocity disper- 
sion and (n) is the mean density. "Larson's laws" were 
derived for a wi de range of str u ctures : from ~ 0.1 pc 
('coherent cores', iGoodman et al.l (|l998l )) up to ~ 100 pc 
(clouds and cloud complexes), and revisit ed many times 
in observational and numerical works ( Solomon et al 



19971: 



19871: iMvers fc Goodman! 19881; Vazquez- Semadeni et al 



Ballcstcros-Parcdcs fc Mac Lo 



iwj 

Dib. Burkert fc Huieiratl 120041 : iHever et al.l |2009h . An 



20021 : 



issue that is still 
exponents a and /3 



not clarified is whether the scaling 
are interdependent - e. g., considering 



comp ressible fluid in a statistical steady state ( Kritsuk et al.l 
l2007f ). or assuming equipartitio n of energies in MC c lump s, 



relations ()Ballesteros- Paredesl 120061 ). 
a recent observational study of giant 



obeying the above 
On th e other hand, 
MCs (|Hever et all I2009D indicates that the coefficient in 
"first Larson's law" scales with their surface density, as not 
expected for purely turbulent obj ects but for ones wherein 
gravi ty plays a significant role (|Ballesteros-Paredes et all 

lion]). 

The structures considered in this work are condensa- 
tions in MCs, generated by turbulent shocks at a given spa- 
tial scale L. Hereafter in this Paper, we label those conden- 
sations clumps, irrespectively of their density n and its ratio 
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to the mean scale densitjQ (n)(L). A turbulent cascade is 
considered within a range of scales 30 > L > 0.1 pc, which 
corresponds to large regions within giant MCs and cloud 
complexes down to dense cloud cores. The upper limit is set 
to provide the molecular phase of the ISM and the applica- 
bility of an assumption of isothermality whereas it is proba- 
bly the very proces s of cloud formation that dri ves the inter- 
nal turbulence (see lKlessen fc Hennebelle]|2010l ). The chosen 
lower spatial limit is a typical transsonic scale at which the 
thermal energy per unit mass becomes comparable to the 
turbulent one. (For further details we refer the reader to 
the implementation of this app roach in the model of the 
stella r initial mass function of IVeltchev. Klessen fc Clarkl 
120111 .) Our study is based essentially o n the results of 
state-of-the-art 3D numerical simulatio ns (|Dib et al 1 l2007l ; 
iKritsuk et "all 120071 ; IShettv et all bOloD . These simulations 
are able to reproduce the inertial range of turbulence only 
for small scales: 0.05 < L < 1.5 pc. To enlarge the scope 
of our study, including clumps generated at scales up to 
~ 30 pc, we re f er as well to the classical 2D simulation 
of iPassot et al.l (|l995l ). with length unit 1 kpc, and its 
implications for the energy balance of clouds and clumps 
|Ballesteros-Paredes fc Vazouez-Semadenilll995h . 

It is natural to expect a mass-density relationship for 
clumps formed in the cloud via turbulent fragmentation. 
This would set a physical link between them and pro- 
tostars. Moreover, a power-law mass-size relationship for 
clouds and clumps in the ISM, term ed "Larson's third law" , 
is confirmed both from observat ions (|Kauffmann et al . 2010; 
lLombardi. Alves fc Ladal |2010| ) and numerical simulations 
(e.g., Shetty et al. 2010). Then, if "Larson's second law" 
(equation[2]) holds for small clumps and cores, although with 
a varying value of the exponent a, a relationship between 
clump mass m and clump density n should be expected. 
The present work is an attempt for theoretical substantia- 
tion of a power-law relationship n oc m x . The exponent x is 
estimated, assuming equipartition between various forms of 
energy within the clumps. 

In Section [5] we present the physical assumptions on 
which our study is based, define the supposed clump mass- 
density relationship and explain the choice of its normaliza- 
tion. Our approach to estimate the power-law exponent x 
as a function of the spatial scale L is described in Section [3] 
In Section [4] we present the results, compare them with ob- 
servations of MCs and discuss the implications of different 
equipartition relations for the clump mass-size relationship. 
Section [5] summarizes the conclusions from this study and 
sketches its possible extensions. 



2 DEFINITION OF THE CLUMP 
MASS-DENSITY RELATIONSHIP 

2.1 Basic physical assumptions and scaling laws 

Our consideration refers to an early stage of the MC 
evolution when turbulence is the key factor for struc- 
ture formation at any spatial scale while gravity takes 
over at small scales. Driven turbulence at a given L is 



1 By use of this ratio, some authors introduce subcategories like 
'voids' and '(dense) structures'. 



fully developed at about one crossing time t cr (L) and 
saturated at ~ 1.5 t cr when a steady state is achieved 
(IPassot et al.l 1 19951; IFederrath. Klessen fc Schmidt! 120081 ; 



iFederrath et al. 2010t ). At the latter epoch, gravity starts 
to play an essent i al rol e in the global evo l ution of the 
cloud (|Elmegreenl |2000| ; IVazquez-Semadenil |2010T ) - its 
total gravitational energy becomes comparable to the 
kinetic and internal energy ( Vazquez-Semadeni et alj [20071 ; 
IZamora-Aviles fc Vazquez-Semadeni! 1201 lh . The fully satu- 
rated turbulence provides a wide inertial range of scales with 
well defined scaling laws for density, velocity and magnetic 
field. Dense structures in the initially diffuse interstellar 
medium form along the magnetic field lines and become 
molecular, self-gravitating and ma gnetically supercritical at 
rough ly the same time sc ale (IVazquez-Semadeni et al.l 
120071 : IVazquez-Semadenil bOloD . with densities of 
a few hundred cm~ 3 or more 
Various simulations show 



2010) 



1995 



1995 



(ICrutcher et al 
(IPassot et al 



.Bal lcstcros- Paredes fc Vazquez-Semadeni 
Clark fc Bonnelll 120051; iHennebelle et al.l 120081 



Baneriee et al.l 120091 ; IShettv et all l2010l ; iDib et all I201C ) 

that such condensations have typical sizes of MCs and 
clumps within them which fall into the inertial range of the 
natal environment. Also, the temporal stationarity within 
the inertial range allows for application of the ergodic 
hypothesis: a time averaged physical quantity can be 
replaced by its ensemble- or spatially averaged counterpart. 

We assume also that the scaling laws do not depend on 
the t urbulent forcing which is a good first- o rder approxima- 
tion (|Federrath. Klessen fc Schmidt! 120081 ; IFederrath et all 
120101 ). The scaling of velocity is generally described by the 
"Larson's first law" (equation [1]) while the d erived value of 
its exponent /3 varies in relatively wide range. iPadoan et all 
l|2006l . I2009T ) found /3 = 0.41 - 0.43 from es timation of the 
turbu lent power spectrum in MCs, while iHever fc Brunt! 
|2004h obtained /3 = 0.65 from observational structure func- 
tions of giant MCs and Monte Carlo modeling. Both works 
combine observational data with synthesized maps but ap- 
ply _diffierent methods. A value /? ~ 0.43 was obtained also 
bv lDib et all(|2007l ). from the correlation of velocity disper- 
sion with size found in their simulations. In this Paper, we 
adopt the alternative values /3 = 0.42 and 0.65, which reflect 
the variation of this exponent and are close to its limit val- 
ues. They are reprodu ced from the numerical simulations of 
IFederrath et all (|201fj| ') by use of the Fourier power spectrum 
and the PCA method, in case of pure solenoidal turbulent 
forcing. 

The density scaling law in the form of "Larson's second 
law" (equation [2]) is derived self-consistently at each scale, 
since it is interdependent with the assumed mass-density 
relationship for MC clumps (see Section [2 . 31 and [3)) . We sup- 
pose that it does not affect the scaling of the velocity. That 
is generally consistent with the very weak c orrelation of the 
local Mach number with density, found by IFederrath et all 
l|2010l l 

The dependence of the magnetic field B on the spatial 
scale is obtained from its relation to the mean density. The 
lat ter is wide l y ado pted to be B oc (n) ' 5 , bas ed on the study 
of ICrutcher! (|l999h . However, as shown in ICrutcher et all 
|2010l l a refined analysis of the available data for MCs yields 
a steeper slope (~ 2/3), in consistency with the theoretical 
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predictions for spherical collapse l|Mestel|[l96r3 ). Therefore 
we use the form: 

0.67 



B = B Q 



flcrit 



(3) 



where n cr it is a critical density at which atomic clouds be- 
come selfjg-avrtati^ and magnetically supercrit- 
ical (|Vazauez-Semadenil l2010l ) . A pronounced scaling law 
B({n)) exists for (n) > n cr i t , where we choose values of 
n C rit = 150 cm -3 (for molecular hydrogen) and Bo = 10 fiG, 
that correspond to it s upper envelope (see Fig. 1 and 5 in 
ICrutcher et "all (|2010h ). 

The clouds and clumps in our consideration are assumed 
to be isothermal with typical temperature T = 10 K. This 
is a good approximatio n for mean densities (n) > 10 2 cm~ 3 



Hennebelle et al1l2008h and corresponds to scales L < 30 pc 
Larso3[l98l|). 

2.2 Clump density distribution and the concept of 
'average clump ensemble' 

The network of interacting shocks generated by su- 
personic turbulent flow causes density fluctuations in 



MCs. As demonstrated from 



many 



numerical sim- 



ulations jKlessenl |2000|; Li. Klessen fc Mac Low! 



Dib et al 



2003 



2008 



Padoan et al. I 120071 ; iKritsuk et al. I 120071 ; _ 
Federrath. Klessen fc Schmiddl200a iFederrath et al.fl20ld ). 
an appropriate statistical description of the distribution of 
the density n per unit volume is a standard lognormal prob- 
ability density function (pdf): 



p(s) ds 



1 



V2ttv 2 



exp 



ds 



(4) 



with standard deviation (stddev) a which depends on the 
Mach number A4 for supersonic flows and determines the 
peak position: 

c 2 

Smax 



~~2 



a 2 = ln(l + b 2 M 2 ) 



(->) 

(6) 

where s = ln(n/(n}) and (n) is the mean density in the con- 
sidered volume which we estimate by use of equation @. 
The quantity b is usually called turbulence forcing parame- 
ter and v aries in the range . 2 — 1.0, depending on the driv- 
ing type llKritsuk et al.ll2007l ; lFederrath. Klessen fc Schmidt] 
l200sl ; iFederrath et alj|2010f i. In the present work we vary 
the values 0.30 ^ b ^ 0.70, depending on the approach and 
the data extrapolation. Such variations of the parameter b 
affect the velocity sca ling slopes in the range 0.42 — 0.48 
l|Federrath et aljfeoiCh . much less than the one covered in 
our study. 

The pdf gives a distribution of density contrasts gener- 
ated through turbulent velocity field. How can one relate a 
density contrast to a specific spatial object with its physical 
parameters (size, mass etc.) like a clump in the ISM? Our 
study is focused not on a set of individual clumps with their 
characteristics but on a group of statistical objects, labeled 
average clump ensemble, representative for all clumps gen- 
erated within a considered volume L s . The average clump 
ensemble at given scale L results from ensemble averaging 
over the variety of Galactic clouds and cloud complexes. 



The peak of the pdf (equation [5]l corresponds to the 
most probable member of this ensemble ('typical clump') 
with density: 



n c = (n) exp(s max ) = (n) exp ( — 



(7) 



Hereafter, we denote the parameters of the 'typical clump' 
by subscript c. The relation (J7J includes dependencies on the 
spatial scale L from equations [1] (through the Mach number 
dependence of the stddev, Eq.[6]) and [2] The range of densi- 
ties in the average clump ensemble is defined as ±cr/2 from 
the most probable density n c . That yields variations in re- 
spect to the mean scale density (n): 



(8) 



(n) exp fs max — < n < (n) exp (s m ax + ^ 

The density of a 'typical clump' differs from (n) within an 
order of magnitude for all considered scales < 30 pc. As will 
be demonstrated in Section [4] its size l c is substantially less 
than the scale L. 



2.3 Clump relationships and their normalization 

As mentioned in Section 1, we assume a power-law clump 
mass-density relationship: 



m 

\n 



n \ I m 
— x In 



\m 



(9) 



where no and mo are units of normalization. Recalling that 
the average clump ensemble and the 'typical clump' are sta- 
tistical objects, equation ([9]) implies a clump size-density 
relationship: 



ct x In 



(10) 



We denote the clump density scaling exponent at a 
given scale by a x with reference to the mean density scaling 
exponent a (equation [2]). Note that these two quantities are 
not necessarily equal - this would be the case if the clumps, 
generated at some scale L, obey the scaling law of the mean 
density. Further in Section 13.21 we assume a self-consistent 
density scaling, requiring a x = a for each L. 

Using equations [9] and 1101 characteristic mass m c and 
size l c are ascribed to the 'typical clump' with density n c : 



m c — mo 



lc = lo ( — - 

\n 



no 



(11) 



(12) 



By use of equations Q , ([9]) and (|10|) , one derives lognor- 
mal mass and size distributions with positions of the maxi- 
mums and widths as follows: 



Smax, m — ^niax/^ , 0~ m — CJ"/|3T 
Smax, 1 = Smax/Oi , O 7 = <j/\a x \ 



(13) 
(14) 



A natural choice of the normalization unit no is the 
mean density at the scale in consideration L: no = (n). 
The choice of the size normalization unit lo should be made 
taking into account typical relative sizes of clumps in respect 
to the scale in consideration. It is appropriate to choose lo 
to be small and proportional to the scale size L: 
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lo = kL 



(15) 



where the dimensionless parameter k could be interpreted 
as mapping resolution of the scale volume. Its appropriate 
value should be chosen to be of order of 10 -2 to achieve a 
distinction of substructures, significantly smaller than the 
spatial scale L and significantly larger than the scale of dis- 
sipation. Further comment on that issue follows in Sect. 14.21 
The relation between the normalization units no, mo 
and lo is obtained from the requirements for volume and 
mass conservation, as shown in the Appendix [X] In view of 
our simple definition of a clump as a condensation of arbi- 
trary density (see Section 1), a scale volume is considered 
to be totally occupied by N clumps. The requirement for 
volume conservation then yields: 



'? ~ Nil exp 



2 (1 
ax- — 



x)(l — 2x) 
2x 2 



(16) 



with the obvious relation between k and N (cf. eauation ll5[) 



2 w {x -l)(l-2x) 
2x 2 



(17) 



On the other hand, the total mass of the scale is M = 
a{p)L 3 = a/i(n)L 3 = apoL 3 where we adopt fi = 2Am p and 
a is a dimensionless parameter of order unity which reflects 
the clumps geometry. The requirement for mass conservation 
leads to: 



M 



Nmo exp ( a 2 x 



2s 2 



(18) 



Then, from equations \W\ and 1181 one obtains the rela- 
tion between the normalization units: 



Polo 

a = exp 

m 



2 1 
ax — 



(19) 



In this work, we assume spherical shapes of clumps and 
spatial scales. Thus a = n/6, but the result is principally 
the same when cubic shapes are considered. 



3 APPROACH TO ESTIMATE THE 
MASS-DENSITY EXPONENT 

The power-law exponent x in equation ([9]) depends on 
the clump physics. It was adopted as a free parame- 
ter in the model of the i nitial stellar mass function of 
IVeltchev. Klessen" fc Clark! I 201lh without theoretical sub- 
stantiation. In this paper we use different approaches to es- 
timate x, based on equipartition relations between the grav- 
itational, kinetic, internal and magnetic energy of the MC 
clumps. 



3.1 Equipartitions of energies 

The energy balance in an object of arbitrary form in the 
ISM is determined by use of the virial theorem (VT). Its 
Eulerian form is given el s ewhere jMcKee fc Zweibel|[l992l : 
iBallesteros-Paredes! 120061 ; iDib et alj|2007f ). In the current 
work we study the effect of relations only between the VT 
volume energy terms since they could be estimated more 
easily from observations. The gravitational, kinetic (turbu- 
lent) energy, thermal (internal) energy and magnetic energy 
per unit volume are defined as: 



W = 

2V 



P 4>dV , E kin = ^ I pu 2 dV, 



8nV 



B 2 dV 



where P is the pressure, <j> is the gravitational potential 
and V is some arbitrary volume. The volume normalization 
is appropriate for clumps since their main parameters are 
derived fSection l2.3[) by using a volume weighted pdf (equa- 
tion [4|. For better readability, we use hereafter simply the 
term 'energy' instead of 'energy per unit volume' when re- 
ferring to the quantities listed above. 

There is no requirement for vir ial equilibrium in our 
treatment. iBallesteros-ParedesI l|2006h showed that such an 
assumption is too strong; it is also not consistent both 
with observations and numerical simu l ations of the ISM 
(|Shadmehri et all 12004 |Dib et al.ll2007l ; IDib fc Kiml 120071 ). 
In fact, MCs seem to be in energy equipartition. In our ap- 
proach we assume relations between W, -Ekin, E t h and -Emag 
of the clumps that hold on statistical terms, for the average 
clump ensemble. These relations are labeled 'equipartitions' 
for simplicity. 

Gravity and turbulence are the main factors for clump 
dynamics at the considered stage of MC evolution (Section 
2.1[) - developed tu rbulence, inherited from the e arlier epoch 
of cloud formation l|Klessen fc Hennebelkj|2010r i. shapes the 
density structure while gravity tends to dominate at small 
scales, in local condensations. Therefore the basic statistical 
equipartition for clumps should be a functional dependence 
like: 



|VK|~/ gk £ kin 



(20) 



where / g k is a coefficient of proportionality. For self- 
gravitating clumps of larger sizes ( clouds), a 'virial- like' rela- 
tion f g k ~ 1 — 2 could be expected (|Vazquez-Semadeni et afl 
12007^ while for clumps of sizes typical for dense cloud cores, 
< 0.1 — 0.2 pc, a more relevant equipartition should include 
the internal energy as well. In order to study different stages 
of clump evolution when gravity takes gradually over, we 
adopt further fiducial values 1 ^ / g k ^ 4. 

Various equipartitions were foun d by 

IBallesteros-Paredes fc Vazquez-Semadem] (| 19951 ) (here- 
after, BV95) from analysis of 2D MHD simulations 
l|Passot et al.|[l995T ). These authors delineated clumps and 
clouds within a large range of sizes 1 < I < 100 pc (see 
their Fig. 2) through isodensity contours at densities 4 to 
16 times larger than the mean density in the simulational 
box (L — 1 kpc). Some relations BV95 discovered are: 

• Equipartitions of the kinetic energy vs. magnetic and 
thermal energy (cf. Fig. 4 in BV95), respectively: 



Ekin ~ E t h 
2-C/kin ^ -E ma g 



(21) 
(22) 



• Equipartition of the gravitational vs. kinetic and mag- 
netic energy: 

\W\ ~2£ kin + £ mag (23) 

In the work of BV95, the energy of gravity due to the ambi- 
ent medium was not taken into account in the term \ W\. We 
allow for moderate gravitational influence from the external 
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cloud in our consideration of the typical clump fSection l3.2l) . 
neglecting possible interaction with the large-scale environ- 
ment. 

The tendency towards equipartition among different 
forms of energy found by BV95 stimulated us to include 
additional equipartitions of the gravitational energy in our 
study: 



\W\ ~ f gt E th 

\W\ ~ /gm-Emag 



(24) 
(25) 



with coefficients of proportionality / gt and / gm . 

The relations (|20|l - (|25p shall be considered as a tool 
to estimate the contribution of differenf forms of energy in 
the clump physics. Whether their use is justified or not is to 
be decided on the base of comparison between our results 
(Section[4]) and the referred observational and numerical ev- 
idence. 



3.2 Clump energies and equipartitions functions 

The equipartition equations from the previous Section can 
be written for clumps generated at a given scale within the 
inertial range, by use of the scaling laws of density, velocity 
and magnetic field strength (Section 12. The energies per 
unit volume are calculated for the 'typical clump' while the 
density range of the 'average clump ensemble' (equation [8]) 
defines their variations. Since the 'typical clump' density 
p c — pn c is proportional (equation [TJl to the mean scale 
density (p)(L) oc L a (a < 0), this approach is comparable 
with the one of BV95 who delineate clumps imposing density 
thresholds, indicative for the corresponding scale (see their 
Fig. 2). 

One gets for the gravitational energy: 



IT*/ 1 3 _ 77lc 



(26) 



where the coefficient z c accounts for the contribution of 
the mass outside the clump to its gravitational e nergy. As 
demonstrated by iBallesteros-Paredes et all (|2009T ) , charac- 
teristic values of z c vary between 1 (vanishing gravitational 
influence) and 2 (strong gravitational influence from the ex- 
ternal cloud). That range is applicable for scales below the 
sizes of giant MC, like in our study. We adopt z c — 1.5 in 
all considered cases. 

The volume term of the internal energy is obtained 
straightforwardly in view of the assumed spherical symme- 
try: 



£th = |-PcT , 
2 p 



(27) 



where 5ft is the gas constant. 

The density of the 'typical clump' p c scales with its 
size according to equation (|12|l . If formation of clumps is a 
turbulent phenomenon, this clump size-density relationship 
should be a self-similar extension of the scaling of density 
according to the "Larson's second law"Q(equation[2]). Hence, 



2 Indeed, iLarsonl lll98lh did not make a sharp distinction be- 
tween clouds and clumps (structures within the clouds) in the 
sample he studied. 



one would expect that the mean density scaling exponent a 
obeys the relation 

3.T 



a = a x = 



(28) 



We derive a and x at given scale L in a self-consistent way as 
follows. We start with some fiducial value a = a^ ' and de- 
rive a solution x = x^ which yields through equation Q280 a 
new value a = or 1 '. Then a new value x = x^ is calculated 
which again is used to calculate a and so forth. The iterative 
procedure converges fast and leads to self-consistent values 
of a and x. 

The magnetic energy could be also expressed in terms 
of p c . In contrast to the treatment of clump density demon- 
strated above, we do not assume a self-similar scaling of B. 
The general magnetic field scaling B(L) is given by equa- 
tion[3] On the other hand, since most of the clumps in a given 
'average ensemble' at scale L are underdense in relation to 
the scale mean density (n) (equation [8} , the magnetic field 
B(l) within the ensembl e is asummed t o obey a scaling law 
of shallower slope 1/2 (|Crutcheu fl999l ). This combination 
of different scaling laws of B provides the maximal possible 
values of -E mag in the corresponding equipartition equations. 
Thus the expression for magnetic energy of a 'typical clump' 
reads: 



B 2 {L) n c 



(29) 



87T (n) 

In fact, such d escription is approp riate for objects with 
higher energies (|Crutcher et alj|20ld ). 

Lastly, to calculate the kinetic energy of the 'typical 
clump', one should implement also the velocity scaling law: 

2/3 

Ekin = TT/OcWc = 7:Pc Uq I t~~ — I i «o=2 km/s (30) 



-Pc Mo I "j 

2 \ 4 pc 



where the chosen normalizatio n units of size a nd velocity 
are based on the simulations of iDib et al.1 (|2007T l. 

Now the exponent x of the clump mass-density relation- 
ship is derived as a zero of the "equipartition functions", 
that correspond to the relations in Section f3. II The way to 



obtain these functions and their form ( equations IB HB6|) are 
given in the Appendix [Bl 



4 RESULTS AND DISCUSSION 

4.1 Mass-density exponent at a given scale 

The exponent x of the clump mass-density relationship de- 
rived at a given scale originates from the density distribution 
of the 'average clump ensemble' (equation [8]) and, hence, 
from the adopted scaling laws of density and velocitjjj. Thus 
x is essentially a parameter of the scale. An observational 
counterpart of this quantity could be derived from mass- 
size relations for (usually) nested structures within MCs, ob- 
tained from extinction maps (e.g. Lombardi, Alves & Lada 
2010; hereafter, LAL10) or column density maps (e.g. Kauff- 
mann et al. 2010): 



M s oc R2 



(31) 



3 Due to the scale dependence of a through the Mach number, 
cf. equation [6] 
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Figure 1. Scale dependence of the mass-density exponent x, derived from cquipartitions of gravitational vs. kinetic energy (lines; 
equation 1201) for the 'typical clump' and adopting alternative values of the velocity scaling exponent. Top: fixed forcing parameter and 
mapping resolution (left: b = 0.50, re = 0.014; right: b = 0.30, ft = 0.01) and varied f g ^ = 1.0, 1.5, 2.0, 4.0 (increasing lincwidth). The 
shaded areas denote the zone of solutions for the average clump ensembles. Bottom: strongly self-gravitating clumps (/gk = 2.0) and 
varying forcing parameter b = 0.30, 0.50, 0.70 (increasing linewidth), simultaneously with 0.03 k ^ 0.001 . A typical solution for a 
scale-dependent re is drawn (grey lines; see text for details). Observational estimates x(-y) according to lLombardi. Alves fc La da (2010) 
for some clouds and cloud complexes (open circles) are plotted for comparison. 



where M s and R s = \JS/ir are the mass and the effective 
radius of a structure with area S. Note that in both men- 
tioned observational studies, the mass-size exponent 7 is 
found to vary with R s . Considering such structures as scales 
of clump formation, i.e. R B ~ L, and taking into account 
the assumed self-consistent density scaling (equation I28[) . 
one derives from equation ()11[) and (|12[) : 



(32) 



7 



Hence one is enabled to compare observational estimates of 
x(-y) (relation I32p with the values of x obtained by use of our 
approach from various equipartitions (equations 1201 - [25]) . at 
different scales L considered as substructure entities of given 
MC. We chose for comparison the work of LAL10 since they 
studied the internal structure of a rich sample of clouds, 
covering a large range of scales (0.01 < R s < 70 pc). The 



latter partially coincides with range of cloud/clump sizes in 
the work of BV95 which inspired our equipartition approach. 

Our method is statistical and, in this sense, universal for 
the variety of Galactic MCs. On the other hand, the results 
depend on several varying parameters: the forcing parameter 
0.30 ^ b ^ 0.70, the velocity scaling exponent 0.42 ^ ft ^ 
0.65 and the mapping resolution parameter n < 0.1. The 
assumption for self-similar density and velocity scaling of 
the average clump ensemble sets a lower limit of the spatial 
scale: the 'typical clump' size is required to fall within the 
inertial range, i.e. its minimal value must be ~ 0.1 pc. Since 
l c is typically about an order of magnitude less than the 
scale of generation for all equipartition methods, relevant 
results could be obtained for L > 1 pc. 

As evident from Fig. Ill II all approaches exhibit a gen- 
eral trend in three aspects: 

• Negative values of x are obtained independently of the 
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Figure 2. The mass-density exponent x, derived from equipartition of the gravitational vs. kinetic and magnetic energy (lines; equa- 
tion I2.il) for the 'typical clump' and adopting alternative values of the velocity scaling exponent. Left: Varied forcing parameter 
b = 0.30,0.50,0.70 (black lines, increasing linewidth), simultaneously with re = 0.035,0.017,0.010, respectively; Right: Fixed forcing 
parameter b = 0.50, varying only the mapping resolution (black lines, 0.01 ^ re ^ 0.016). Other designations are the same like in Fig. [l] 



approach. Their range at all considered scales is in agree- 
ment with the estimates from the most recent observa- 
tions (Kauffmann et al. 2010, LAL10). Only the method 
Skin ~ £"th yields positive x as well, clearly outside the zone 
of observational data (Fig. [3]). 

• Negative values of x can vary from —0.2 down to —2.5 
while the specific functional behavior (convexity/concavity, 
steepness) depends on the chosen set of values of the param- 
eters b, P and k. When (6, fS) are fixed, increase of n leads 
to a shift of the curve x(L) downwards (i.e. decrease of x), 
retaining its shape. On the other hand, increase of b or /3 
steepens the curve at larger scales. Positive x (obtained in 
the case E^ n ~ E t h) are virtually independent on L and 6; 
the curves x(L) are slightly shifted vertically by variations 
k and (3. Choice of a scale-dependent k generates shallow 
curves x(L), in agreement with the structure of some larger 
cloud complexes. 

• The coefficients of proportionality in equations 1201 1241 
and 1251 affect the curves x(L) in different ways: variations 
of / g k lead to vertical shift; the steepness of the curves is 
highly sensitive on variations of / gt ; while variations of / gm 
do not influence significantly the functional behaviour. 

We sorted out the physically meaningful solutions for x 
by two criteria: i) comparison with the observational data 
2(7) , and ii) comparison with the expected range of the basic 
clump parameters from simulations like size, mass, density 
and mean density of the scale. Both criteria were applied 
simultaneously. 

The results obtained from the equipartition of the 
clump gravitational vs. kinetic (turbulent) energy fit well 
with some of the LAL10 data (Fig. [TJ, independent on the 
chosen velocity scaling. Variation of the factor / g k (top pan- 
els), which may correspond to different evolutionary stages 
of a clump, allows for fitting the observed structure of some 
large cloud complexes like California, Orion A and p Ophi- 
uchi, while it is problematic to attain agreement with a 



'shallow structure' in terms of x (i.e. 7 ~ const; Pipe neb- 
ula) or with a too 'steep structure' (Lupus 1). The results 
are similar for strongly self-gravitating clumps (/ g k = 2.0, 
Fig- HI bottom panels) when the turbulent forcing parameter 
is varied. The effect from variation of the mapping resolu- 
tion parameter k (not illustrated) is a moderate steepening 
of the curves x(L) at larger scales without essential change 
of their shape. Some additional opportunity to fit shallowly 
structured clouds like Pipe nebula gives the choice of a scale- 
dependent mapping resolution parameter (Fig. [T] bottom) 
k oc L~ a , where 0.4 < a < 1.3 for all equipartitions. That 
is a reasonable approach since the wide density and mass 
distributions produced at large scales (cf. equations [1] and 
[6]l require a better mapping resolution (decrease of k) in or- 
der to trace the substructures. However, such an approach 
needs refinement, taking into account the diversity of clump 
shapes. 

In view of the range of possible values of / g k and b, gen- 
eral agreement between the cloud structure, predicted from 
the equipartitions \W\ ~ /gk-Eidn, and the LAL10 data for 
some clouds is not enough to determine the energy balance of 
the 'typical clumps' in them. Rather, we aim to demonstrate 
the capacity of our statistical approach, combined with the 
equipartition equation to derive x. To model the cloud struc- 
ture in more detail, one needs additional data about: i) the 
evolutionary status of the cloud which restricts the applica- 
bility of our concept (Sect. [2~Tjl ; ii) the physical conditions 
in the cloud that determine and/or could be implied from 
its clump mass function. The latter is to be demonstrated 
in a subsequent paper. 

The addition of a magnetic term in the equipartition 
of the gravitational vs. kinetic energy (equation I23[) for 
strongly self-gravitating clumps (/ g k = 2) and k = const 
generates curves x = x(L) with a different shape which can- 
not fit the observational data (Fig. [2}. On the other hand, 
adopting a scale-dependent k as mentioned above yields so- 
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Figure 3. The mass-density exponent x, derived for equipar- 
titions of the kinetic energy (lines; equations I21I and I22[l . Top: 
Ekin ~ -Eth : effects of constant (b = 0.60, ft = 0.024, red line) 
vs. scale-dependent mapping resolution (b = 0.33, grey line) are 
juxtaposed. Note the obtained positive solutions x pa const > 
(b = 0.33, k = 0.1); Bottom: 2E^- ln ~ E mag : only solutions with 
a scale-dependent mapping resolution are possible (k oc L _1 ). 
The effect of varying the forcing parameter b = 0.30, 0.50 arc 
illustrated (increasing linewidth). 



lutions that are consistent with the structure of the large 
cloud complexes; in contrast to the result when a magnetic 
term is missing (Fig. [1] bottom). That points to the impor- 
tance of clump geometry when the magnetic field is consid- 
erecQ. 

Equipartitions of the kinetic energy (equations 1211 and 
I22p yield solutions with anomalous behavior of x; i.e. either 
lead to increasingly steeper slope or to extremely small val- 
ues at small scales (Fig. [3] top). Equipartition Skin ~ E t h 
produces a family of curves with positive x, practically inde- 
pendent on scale. However, those solutions fail to meet the 
criteria we adopted: they lie outside the observational range 
and correspond to an unrealistic set of clump parameters, as 
will be shown below. Negative x that lie in the domain of the 
LAL10 data are obtained in a very limited range of spatial 
scales 5 < L < 10 pc, but are clearly unrealistic. The results 



Note, however, that in the assumed evolutionary time frame 
of our study clump geometry is determined mainly by gravity. 



for the case 2Emh ~ E ma g are similar in terms of appli- 
cation (Fig. [3] bottom). Evidently, the clump mass-density 
relationship cannot be explained only by use of equiparti- 
tions of the kinetic energy, without invoking gravity. That is 
consistent both with time frame of our study: the transition 
between the turbulent epoch and the stage of hierarchical 
gravitational collapse, and with the supposed spherical sym- 
metry of the clumps. Turbulence is the main factor for the 
formation of clumps but gravity is crucial to determine their 
fundamental physical characteristics. The latter suggestion 
is supported by the results obtained from equipartitions of 
the gravitational vs. thermal and vs. magnetic energy (equa- 
tions [24] and 1251 Fig. [4| . The agreement with observations 
of large cloud complexes (California, p Oph) is remarkable 
for choice of a scale-dependent k oc L~ a , except in the case 
\W\ ~ E mag, P — 0.65. The choice of constant k is able 
to explain 'steeply-structured' clouds like Lupus 1, only by 
use of the classical Jeans approach, i.e. gravity vs. internal 
energy. The solutions for / gt = 2 in equation [24] are highly 
sensitive to the choice of j3 and k, retaining the shape of the 
curves, while those for / gm = 2 in equation [2S] practically 
coincide with those for / gm = 1. (Therefore these cases are 
not illustrated in Fig. [4]) 

From the point of view of the current understanding of 
MC physics, we suppose that the better agreement of the re- 
sults for the equipartitions of W with observations of LAL10 
could be explained with the evolut ionary status of the clouds 
in wh ich clumps are generated. iBallesteros-Paredes et all 
(|201lh suggest that the MCs undergo two stages: an epoch of 
hydrodynamic turbulence, followed by an epoch of a hierar- 
chical gravitational collapse. The clumps considered in this 
work are typical objects existing in the transitional epoch 
between them. 



4.2 Mass-density exponent from global mass-size 
relationships 

Clumps observed in MCs are generated by turbulent flows 
at different scales. Nevertheless, irrespectively of the scales 
individual clumps originated from, their masses and sizes 
correlate strongly, as confirmed both by observations and 
numerical simulations (|Larsonlll9STl ; iKauffmann et al.ll2010l ; 
IShettv et a"i]|2010h . That is an indication of a global mass- 
size relationship, with a power-law exponent 7 g iob- From ob- 
servational or numerical estimates of 7 g iob , one can derive by 
use of equation (|32|1 a global mass-density exponent a;(7 g i b). 
In our treatment, such an approach is realized in two steps: 
a) composition of size-mass diagrams for the set of 'typi- 
cal clumps' for all scales L and calculation of 7 g iob from a 
power-law fit of the size-mass correlation; b) derivation of 
a;(7 g i b). The latter quantity is representative for the whole 
SF process at scales 0.1 < L < 30 pc and thus can be a clue 
to derivation of the clump mass function. Note that even if 
7 g iob is found to vary with L (i.e. no single mass-size power- 
law relationship is applicable), a;(7 g iob) is not expected to 
exhibit the same functional behavior like x. 

Composite size-mass diagrams, obtained from the 
equipartition relations (equations 2012511. are plotted in Fig. 
15151 together with the results of IShettv et al.l (|2010l ) from 



simulations. The result for the equipartition of the gravita- 
tional vs. kinetic energy is illustrated in Fig. [5] In the left 
panel, we demonstrate the effect of varying the factor / g k. 
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Figure 4. The mass-density exponent x, derived from equipartition of the gravitational vs. internal (equation I24t and magnetic energy 
(equation I25|l . Top (\W\ ~ Ethi fgt = 1): Solutions for b = 0.50,0.70 (increasing linewidth) for constant (black lines; n = 0.023,0.010 
/left/, k = 0.013,0.007 /right/, respectively) and scale-dependent mapping resolution parameter (grey lines) are plotted; Bottom (\W\ ~ 
£ma g , fgm = 1): solutions for b = 0.30,0.70 (increasing linewidth) for constant n (black lines; k = 0.030,0.008 /left/, k = 0.017,0.006 
/right/, respectively). 



The larger it is, the more clump mass is required to ful- 
fil the equipartition equation. Applying the classical Jeans 
criterion under the assumption of spherical clump shape, 
one sees that part of the clumps are gravitationally unsta- 
ble (filled symbols). Most of them are generated on larger 
scales. There is a lower scale limit for producing u n stable 
clumps, see e.g. Fig. 1 in lVeltchev, Klessen fc Clarkl (|201ll ) 
and the comments on it. Remarkably, the points for all 'typi- 
cal clumps' and even average clump ensembles (the latter are 
not plotted for clarity) are aligned without dispersion. This 
result is not trivial since x, which determines the size-mass 
relationship for the 'typical clump' in an individual ensem- 
ble, is a complex function of a, j3, b and k which includes 
exponential and power-law terms (cf. equation ll9p : 

m c = Ipoll L 8 -/' 1 -' exp ( - (33) 

The slope of the obtained size-mass relationship de- 
pends exclusively on the choice of the velocity scaling pa- 
rameter and is insensitive to / g k, the type of turbulent forc- 



ing and the mapping resolution. This is demonstrated in 
the right panel of Fig. in the case of strongly gravitating 
clumps as the parameters b and n are varied simultaneously. 
Results for a scale-dependent k exhibit similar slope but es- 
sentially lower sizes and masses. 

As suggested from numerous observations of dense 
structures in MCs, they exhibit power-l aw mass-size rela- 
tionships of slope abo ut 2 and larger (e.g.. lOnishi et al ll99d : 
iTachihara et alj|2002t ). Our 'typical clumps' are expected to 
follow similar behaviour as they span the same range of mass 
and size. As seen from Fig. [5] this is indeed the case with 
the equipartition |W| ~ /gk-Ekin- We find 7 g iob ~ 2 for the 
whole range of values of the velocity scaling parameter /?. 
This result could be interpreted in terms of 'typical clump' 
gravitational potential <f> c , since in our approach for deriva- 
tion of x all equipartition functions (see Appendix [Bj that 
include gravity can be expressed through <p c (i.e. include 
density to the first power). If a global size-mass relationship 
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Figure 5. Global size-mass relationships, obtained from the equipartitions of gravitational vs. kinetic energy for /3 = 0.42 (red symbols), 
f) = 0.65 (blue symbols) and adopting a scale-free ft. Left: variyng factor / g it = 1.0, 1.5, 2.0, 4.0 (increasing pointsize) and sets of (ft, k) 
like in the top panels of Fig. [T] Right: strongly self-gravitating clumps (/ g k = 2.0) and sets of (b, k) like in the bottom panels of 
Fig. Q] Gravitation ally unstable clumps are denoted by filled symbols. The zone, populated by the PPP clumps from the simulation of 
Shct tv et al.l i2010f ). is shown (shaded area) for comparison. 



m c oc l2 elab exists, the dependence of the potential on the 
clump size must be 

(f> c oc rric/lc oc l2 slab 1 . 

In view of the abovementioned studies of MCs and cloud 
complexes, wherein 7 g iob > 2 was found, <j> c of the clumps 
should grow linearly or steeper with their size. On the other 
hand, 7 g iob = 1 is the lower limit of observationally verifiable 
clump mass-size relationships. 

Inclusion of a magnetic term in the equipartition by use 
of equation 1231 yields significantly a shallower slope (Fig. [5]). 
It points to the critical role of the magnetic energy in the 
clump energy balance. This can be seen also in Fig. [8] (bot- 
tom) where the results from the equipartition of the gravita- 
tional vs. magnetic energy are shown. Proper description of 
clump geometry is crucial in all cases with a magnetic term, 
as expected in view of the fact that the magnetic field is 
not spherically symmetric. Solutions with a scale-dependent 
mapping resolution (k oc L~ a ) yield strange clump families 
and fail to reproduce the expected size-mass relationship 
(Fig. [6] and [U bottom). 

The gravitional instability of the 'typical clumps' is es- 
timated by use of the Jeans criterion, comparing m c and 
mj(fic). Excluding the ratio n c /no in equations II 11121 and 
replacing mo from equation 1 191 one obtains: 

n i3 ( 2 x ~ 1 \ /„ ,i 

m c = -ttpJc exp a x (34) 

6 V x J 

This relation holds for all members of the 'average 
clump ensemble'. Obviously, m c depends on 'typical clump' 
size and density in a more complex way than naively ex- 
pected, due to the statistical nature of the studied objects. 
Thus the Jeans criterion for instability is different for each 
chosen equipartition. 

The equipartitions of the kinetic energy (equations [21] 
and 122 p do not yield a pronounced mass-size relationship 



(Fig. [7J. The equipartition E^ n ~ E t ^ yields negative so- 
lutions x within restricted spatial scale range (cf. (Fig. [3} 
and predicts clumps of very low sizes. On the other hand, 
positive solutions for x ~ 0.25 — 0.5 are not supported by 
LAL10 data (Fig. [3J) and correspond in general to small 
clumps of very low mass. However, the upper envelope of 
this clump population in the m — I diagram apparently falls 
within the zone of the so called 'coherent cores' (0.04—0.1 pc, 
~ 1 Mo). Indeed, this picture is to be expec ted from the ob- 
serva tional definition of the coherent cores (|Goodman et al.l 
1 1998h as objects wherein the contribution of the thermal ve- 
locity becomes comparable to the turbulent one, Uturb ~ ^th- 
The latter relation corresponds exactly to the equipartition 
Skin ~ E t h (cf. also equation IB3|I . Interestingly, the only 
physically meaningful solutions for the equipartition in ques- 
tion are obtained for the coherent cores. 

The clump families, obtained from the equipartition 
2-Ekm ~ Bma g exhibit anomalously steep slopes 7 g iob that 
lead to values of x > 0.75. As evident from Fig. such 
families could explain clumps in restricted ranges of size 
(I ~ 0.6—1 pc) and mass. Those clumps are generated within 
restricted ranges of scales as well: L > 7 pc. Equipartitions 
of the clump kinetic energy, without accounting for gravity, 
appear insufficient for description of clumps in broad ranges 
of mass and size. 

As already pointed out, the case of equipartition of the 
gravitational vs. the turbulent (kinetic) energy - whether 
in strongly gravitating clumps or not - fits well both 
with observat i ons ( LAL10, Fig. [TJ) and with simulations 
l|Shettv et alj (|2010f ). Fig. |3J. That hints at the substan- 
tially lower weight of the thermal and magnetic energy 
in equipartitions between volume energy terms. Mass-size 
relationships, derived from equipartitions of _E t h or i? mag 
with \W\, are instructive in that aspect (Fig. [HJ. The 
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Figure 6. Global size-mass relationships, obtained for cquiparti- 
tion of the gravitational vs. kinetic and magnetic energy (equa- 
tion !23l l. The designations are the same like in Fig. [5] The effect of 
varying turbulent forcing parameter 0.30 b ^ 0.70 (increasing 
pointsize, /3 = 0.42) and of adopting a scale-dependent re (grey 
symbols, /3 = 0.65) are shown. 
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Figure 7. Global size-mass relationships, obtained from the 
equipartition of the kinetic vs. internal energy (squares, equa- 
tion [2TJ and vs. magnetic energy (circles, equation [22} . In the 
former case, only clumps that correspond to positive x are shown. 
The sets of b and re are the same like in Fig. [3] (top and bottom). 
Other designations are the same like in the previous figures. 



equipartition \ W\ ~ /gt-Eth (equation I24|) recovers perfectly 
the classical Jeans approach, i. e. the 'Jeans-like case' in 
IVeltchev. Klessen fc Clarkl (j201lf l. Under the assumption of 
isothermality, mj oc n~ 1//2 oc m m oc mj, than 

m oc /, i.e. 7 g iob = 1 (Fig. [HI top). In other words, the case 
\W\ ~ /gt-Eth produces clumps of mass, proportional to the 
ensemble averaged Jeans mass. Appropriate choice of the 
factor fgt > 2 could determine the domain of gravitational 
instability. The equipartition \ W\ ~ / g m£mag yields a signif- 
icantly shallower slope of the mass-size relationship, similar 
to the case when i? mag and E^- ln are both taken into consid- 
eration (Fig. [6}. It confirms the critical role of the magnetic 
energy in the balance between clump volume energy terms. 



We point out again the significance of the clump mapping 
in this case, the choice k oc L~ a leads to unrealistic results. 
We attribute them to the specific clump geometry shaped 
by the presence of magnetic fields. 



5 SUMMARY AND CONCLUSIONS 

This work is an attempt for theoretical substantiation of the 
clump mass-density relationship n oc m x and for a plausible 
estimation of its exponent x. Due to the turbulence and the 
dynamics of the ISM, such a relationship has a statistical 
significance and therefore a statistical approach was chosen. 
The time frame of our consideration refers to the transition 
between two main stages of the MC evolution: the predom- 
inantly turbulent epoch and the period of subsequent hi- 
erarc hical gravitational collapse (|Ballesteros-Paredes et aLl 
2011). At that timescale, condensations in the initially dif- 
fuse, atomic ISM become molec ular, gravitationally boun d 
and magnetically supercritical l|Vazquez-Semadenil |20T0T ). 
The turbulence in the newly formed MCs becomes essen- 
tially supersonic, due to the dra s tic decrease of tempe rature 
(|Vazquez-Semadeni et"al] 120071 ; iBaneriee et"alll2009| ). The 
turbulent inertial range covers spatial scales from several 
tenths to several decades of parsec. We assume steady state 
and validity of the ergodic hypothesis. Using the scaling laws 
for velocity, density and magnetic field and equipartition re- 
lations for various forms of energy for a statistically averaged 
('typical') clump at different scales L, we obtain equations 
to derive the mass-density exponent x. 

The obtained results could be summarized as follows. 

1. Generally, the solutions x correspond, both as range 
of values and as functional behavior, to the (sub)structure of 
MCs, recently studied by LAL10 by use of extinction maps. 
However, this agreement depends substantially on the cho- 
sen equipartition relation. The general conclusion is twofold: 
i) equipartitions which do not include gravitational energy 
are not relevant for derivation of x, and, ii) in contrast to 
gravitational and kinetic energy, magnetic and thermal en- 
ergy do not play significant role for the clump mass-density 
relationship. 

Good fits to the observational data of LAL10 are ob- 
tained by variations of the turbulent forcing parameter b 
and the mapping resolution k and are not sensitive to the 
choice of velocity scaling law (the value of j3). That sug- 
gests a connection between those parameters and the spe- 
cific structure and evolution of the clouds and, possibly, to 
an universally applicable range of variations of /3. However, 
to model the cloud structure in more details, one needs ad- 
ditional data concerning the evolutionary status of the cloud 
and physical conditions therein. Our aim in this work is sim- 
ply to demonstrate the capacity of the proposed statistical 
approach, combined with equipartition equations to derive 
x. Finally we point out that variations of the clump param- 
eters within statistically fiducial ranges ('average clump en- 
semble') also cover the diversity of MC structure, exhibited 
by the observational analysis of LAL10. 

2. The composite mass-size diagrams, constructed for all 
'typical clumps', generated at different scales, yield a global 
value of mass-size exponent 7 g iob and hence a global mass- 
density exponent a:(7gi b). The latter quantities are repre- 
sentative for the diversity of Galactic MCs and of their sub- 
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Figure 8. Global size-mass relationships, obtained from the equipartition of the gravitational vs. internal (top) and magnetic energy 
(bottom). Squares denote the case / g t = 1 (or, / gm = 1) and circles - the case / g t = 2 (or, / gm = 2). Symbols with black borders are 
examples for solutions k oc L~ a . 



structures. An immediate conclusion from the clump mass- 
size diagrams is that an equipartition relation between the 
clump gravitational and kinetic energy provides realistic es- 
timates both of x and of x("f s iob)- This lends additional sup- 
port to our previous conclusion about the secondary im- 
portance of the clump magnetic and thermal energy for the 
clump mass-density relationship. On the other hand, the im- 
plied geometry in the description of magnetic field is crucial 
when a magnetic term is present in the equipartition rela- 
tions. 

3. In contrast to the local value x, the behaviour of the 
global mass-density exponent x("/ g i b), derived from differ- 
ent equipartitions which include gravitational and kinetic 
energy (Fig. [5J, does not depend on the turbulent forcing, 
on the mapping resolution and on the equipartition coeffi- 
cient. It leads us to conclude that the obtained estimates 
of the global slopes of the mass-size and mass-density rela- 
tionships are universal for the variety of Galactic MCs and 
their substructures. However, those slopes depend notably 
on the choice of the velocity scaling law. This finding could 
be interpreted in terms of the prehistory of MC formation. 



Turbulence in the MC is inherited from the natal diffuse ISM 
and hence its power spectrum is preliminarily determined. 
We suppose that the prehistory of MC formation affects the 
clump mass function as well, through the value of 7 g iot>- 

4. The analysis of our results demonstrates the impor- 
tance of the different forms of energy in the dynamics of gas 
clumps on the verge between the purely turbulent epoch in 
the MC evolution and the stage of hierarchical gravitational 
collapse. The key factors for clump formation and struc- 
ture are gravity and turbulence. Our assumption of spheri- 
cal symmetry of the clumps implicitly diminishes the signif- 
icance of the magnetic field. The role of the thermal energy 
would increase if spatial scaling of the temperature is in- 
cluded in the description. 

Possible future extension of this work is to implement 
spatial scaling of the gas temperature and to improve the 
geometrical description of clumps, retaining our basic phys- 
ical assumptions. The estimates of the global clump mass- 
size slope could be used to develop a model of the clump 
mass function as a key to the stellar initial mass function. 
It would be also interesting to extend the time frame of our 



Mass- density relationship in molecular cloud clumps 13 



consideration, including earlier (non-gravitational) and later 
(gravitationally dominated) stages of the MC evolution. 
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APPENDIX A: RELATIONS BETWEEN THE 
NORMALIZATION UNITS 

If Ni is the number of clumps of size I and Si = \n(l/lo), the 
requirement of volume conservation gives: 



v = ajZN i = N a i 3 ( N i/ N ) = N al3 p( l ) 



(Al) 



where p(l) = p(si) dsi is the probability (fraction) of clumps 
of size I and TV is the total number of clumps, generated at 
L. Thus, since p(si) dsi = p(s) ds (cf. equations [13} 14), one 
obtains 



V 



aN J l 3 p( Sl ) da, = aNl 3 J p(s) ds 



Using equation 1101 one gets: 

r, Alii 

V 



aNlt 



— exp 

n 



1 / S Sma 

2 



ds 



Further, from n/no — exp(s) and by rewriting the integral 
for variable y = (s — s m ax)/V2cr 2 , we obtain: 



v aNlt fl — x 
V = exp ( 



exp(cj/ -y )dy , 
1-x 



c = V2a 2 



(A2) 



The right-hand side of this equation contains a Gauss inte- 
gral with an analytical solution. Taking into account that 
V = aL 3 , we derive finally equation [16] 

The requirement for mass conservation yields equation 
1181 in analogical way, starting from: 



— p{s) ds (A3) 



M ~ N / mp(s m ) ds m = Nmo 



APPENDIX B: EQUIPARTITION FUNCTIONS 

These functions are used to derive the clump mass-density 
exponent i at a given scale. They are obtained from the 
equipartition relations, given in Section [3.1l through the fol- 
lowing procedure: 

- Expressions for different energies (|26l - 130|) are substi- 
tuted in the equipartition relations (j 201) - (|25|) 

- The quantities l c and m c are excluded from the equa- 
tions via substitution with p c by use of equations (|11[) and 

- The 'typical clump' density p c is expressed through 
equation (0. 

- Clump mass normalization unit mo is expressed 
through po and lo fequation ll9|) . 



Bl Gravitational vs. kinetic energy 
From equation (1201) : 

'Ax + 2 ar-1 



Q g k(a;) 



— ZcGpoh exp 
5 



6,r 



B2 Gravitational vs. kinetic and magnetic energy 
From equation (|23l) : 

'Ax + 2 x-1 



Qgkmag(^) 



t ,a 
— ZcGpolo exp 
5 



~ W ° ( I^J 6XP 



(i.r 



2 2/3 + (3 - 2/3)a; 

(7 X - 



(i.r 



87rpo 



exp 



(T 

~2 



(B2) 



B3 Kinetic vs. thermal or magnetic energy 

From equations ([21) and ([22)) : 

2 ^ 2/3 + (3 - 2P)x 
6x 



'kth(x) = "0(7-^-) exp ( - a > 



A pc 



■ 3exp( ~tJ — 



2;i 



?kma g (a;) = UqI ^-77 ] exp ( — a' x 



4 pcj 



5" / of 
' 87rpo 6XP V 2 



(B3) 

2 2/3 + (3 - 2/3)a; 



6.r 



(B4) 



B4 Gravitational vs. thermal or magnetic energy 
From equations ([24]) and ([25]) : 

'Ax + 2 x-1 



Qgth(a;) 



^ n 72 
— z c &polo exp 

5 



3, / a* 
--/ st exp( - 



6x 



(B5) 



Qgmag(x') = — 2: c Gp0^exp 



2 / Ax + 2 x-1 
a x 



(:>.<: 



B 

1 87rpo 



exp 



a 
~2 



(B6) 



The obtained equipartition functions are listed below. 



